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ABSTRACT 

We report new correlations between ratios of band intensities of the 15-20 pm emission bands of 
polycyclic aromatic hydrocarbons (PAHs) in a sample of fifty-seven sources observed with Spitzer/IRS. 

This sample includes Large Magellanic Cloud point sources from the SAGE-Spec survey, nearby 
galaxies from the SINGS survey, two Galactic ISM cirrus sources and the spectral maps of the Galactic 
reflection nebulae NGG 2023 and NGG 7023. We find that the 16.4, 17.4 and 17.8 pm band intensities 
are inter-correlated in all environments. In NGG 2023 and NGG 7023 these bands also correlate with 
the II.O and 12.7 pm band intensities. The 15.8 pm band correlates only with the 15-18 pm plateau 
and the 11.2 pm emission. We examine the spatial morphology of these bands and introduce radial 
cuts. We find that these bands can be spatially organized into three sets: the 12.7, 16.4 and 17.8 
pm bands; the 11.2, 15.8 pm bands and the 15-18 pm plateau; and the 11.0 and 17.4 pm bands. We 
also find that the spatial distribution of the 12.7, 16.4 and 17.8 pm bands can be reconstructed by 
averaging the spatial distributions of the cationic 11.0 pm and neutral 11.2 pm bands. We conclude 
that the 17.4 pm band is dominated by cations, the 15.8 pm band by neutral species, and the 12.7, 

16.4 and 17.8 pm bands by a combination of the two. These results highlight the importance of PAH 
ionization for spatially differentiating sub-populations by their 15-20 pm emission variability. 

Subject headings: astrochemistry, infrared: ISM, ISM: lines and bands, ISM: molecules, molecular 
data, techniques: spectroscopic 


1. INTRODUCTION 


Polycyclic aromatic hydrocarbons (PAHs) are t he most 
comm on polyatomic molecules in the universe (Tielens 
2013) and are thought to produce prominent infrared 


emission bands at, e.g., 3.3, 6.2, 7.7, 8.6, II.2 and 12.7 
pm. These emission bands are seen in many differ¬ 
ent environments, including Hll regions, reflection neb¬ 
ulae (RNe), planetary nebulae (PNe), and the diffuse 
interstellar medium (ISM). They are attributed to in¬ 
frared (IR)-active vibrational modes, generally stretching 
or bending modes of G-G and/or C-H bonds. Bending 
modes of the carbon skeleton, denoted as C-G-G modes, 
produce weaker but identifiable emission features at 15.8, 
16.4, 17.4, and 17.8 pm, typically perched atop a broad 
17 pm plateau (e.g., Monton et al. 2000 Van Kerckhoven 


et al.||2b00[|Van Kerckhoven||2002 ). 

PAH emission bands can be separated into three cate¬ 
gories. The 3-15 pm bands arise from nearest-neighbour 
vibrations, such as the 3.3 pm C-H stretch or the 8.6 
pm in-plane C-H bend. These nearest-neighbour vi¬ 
brations are common to all PAHs, and their collec¬ 
tive emission features at similar wavelengths leads to 
blended emission. As a result, the emission intensities 
are high and they are readily observed in astronomi¬ 
cal environments. In contrast, vibrations of the entire 
molecule (typically occuring at > 20 pm, such as “drum- 
head” modes) are extremely molecule-dependent (e.g 


Van Kerckhoven et al.||20Q0| [Boersma et alT 2010| Ricca 

et al.|20IQf Boersma et al.|20Iip . The intensities of these 
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bands are very low and they are difficult to detect; in 
fact, none have, to date, been identified in space. The 
15-20 pm bands however represent an intermediate be¬ 
tween the strong 3-15 pm bands and the rich but weak 
20-f pm bands; they arise from vi brational modes involv¬ 
ing a subset of the inolecnle (e.g.. Van Kerckhoven et al 


a. 'i’hey' 
50 - 200 


2000 Peeters et al. 2004 Boersma et al.'2010| 
are generally t hought to arise in PA Hs witn 
carbon atoms (Boersma et al.||2010). The plateaus upon 
which the PAH bands sit are thougnt to ari se from emis¬ 
sion by large PAHs or P AH clusters (e. 
et al.||1989 Tielens||2008 ). 


Allamandola 


The 3-15 pm bands are known to have a series of corre¬ 
lations between the intensities of specific bands, typically 
attributed to common properties (primarily ionization 
state). The 15-20 pm PAH bands however have thus far 
not exhibited any correlations amongst themselves (e.g.. 
Smith et al. 20Q7b| Boersma et al. |2010| Ricca et al. 
2010[ [Peeters et arf|2012[ hereafter referred to 

I), though there are results showing the 16.4 


as paper 
pm band 


corre lates with the 6.2 and 12.7 pm bands ([Boersma et al 
2010 paper I). Previous studies have had limited sample 
sizes of moderate quality or focused on single objects. 
Here we investigate correlations in a large sample of ob¬ 
jects to examine the behaviour of the 15-20 pm bands, 
and their relationship to the 3-15 pm emission features. 
The paper is organized as follows: the observations and 
data reduction are presented in Section 2. The methods 
used to analyze the spectra are summarized in Section 
3. Results, including correlation plots and spatial maps, 
are presented in Section 4. We interpret these results in 
Section 5 and summarize our findings in Section 6. 


2. OBSERVATIONS AND DATA REDUCTION 
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2.1. Target selection 

We obtained spectra of a variety of objects from the 
SAGE -Spectroscopy (SAG E-Spec) Spitzer legacy pro¬ 
gram (Kemper et al. and sources from the Spitzer 

Infrared JN earby Galaxies Survey (SINGS, |Kennicutt 
et al.|2003|) . We chose sources displaying clear FAH emis¬ 
sion in the 15-20 pm region. Such emission can be identi¬ 
fied by a wide 17 pm complex (or plateau) and narrower 
features at 15.8, 16.4, 17.4 and 17.8 pm. Sources more 
extended than 13.5" in the SAGE-Spec sample were ex¬ 
cluded, as this corresponds to the approximate size of the 
SL aperture at 11.2 pm (which is used for normalization 
when comparing band intensity ratios). Only sources 
with a 3(7 detection of their 15-20 pm emission compo¬ 
nents (or a subset) are included. We also applied this 
measurement criterion to the “MIR” (5-15 pm) PAH 
bands. Eifty-seven sources were identified as meeting 
these selection criteria, including two ISM cirrus sources. 
The sources in our sample are presented in Table We 
also analyze a high-resolution spectral map of the RN 
NGG 7023 (Eigurepl previously analyzed by Rosenber; 


2014a 


et al.|20TT] Berne fe^Tielens 2Q12[ Boersma et al. 


2015). We compare our results to that of NGG 
I and Peeters et al. 


2015, in prep., here¬ 
after refe rred to as paper II). The 15-20 pm measure¬ 
ments of Boersma et al. (2010), comprising 15 sources, 
are also included for comparison. Our total sample in¬ 
cludes carbon-rich evolved stars, young stellar objects 
and Hll regions in the Large Magellanic Gloud, nearby 
galaxies. Galactic RNe, and other ISM sources. 

2.2. Observations 

The spectra we analyze were acquired wi th the In¬ 
frared Spectrograph (IRS; H ouck et al. 20041) on board 
the Spitzer Space Telescope (|Werner et al.||2004|). The 
IRS instrument provides high- and low-resc)lution spec¬ 
troscopy in the mid-IR regime (5 — 40 pm). Low- 
resolution spectra (R ^ 60 — 130) were obtained with 
the Short-Low (SL) and Long-Low (LL) modules, which 
provide coverage from 5 — 14 pm and 15 — 40 pm, respec¬ 
tively. We also analyze a high-resolution {R ^ 600) spec¬ 
tral map of NGG 7023 and 2 ISM objects acquired with 
the IRS/Short-High (SH) module, covering the wave¬ 
length range 9.9 — 19.5 pm. These SH data were ob¬ 
tained from the NASA/IPAG Spitzer Heritage Archiv^ 
(AORkey: 3871232, PI: Giovanni Fazio). Previous stucK 
ies that have examined these data are t hose of|Sellgr en 


let al. 


J2Q07|) ; Rosenberg et al. ( |2011| ); [B 




boersma 


et al 


2.3. Reduction 

The SAGE-Spec data included in this study were ob¬ 
tained in reduced form. The project an d the data re- 
( luctM n process are described in detail by |Kemper et al.| 
(2010); we briefly review the reduction steps here. The 
Spitzer Science GenteiH (SSG) data-reduction pipeline 
(version S18.7) was us^ to produce flat-field corrected 
images. Background subtraction was completed using 
on- and off-source nod positions. Rogue pixels were re¬ 
moved and replaced using the imclean algorithm, part of 

^ http://sha.ipac.caltech.edu/applications/Spitzer/SHA/ 

^ http://ssc.spitzer.caltech.edu 
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Figure 1. Hubble Space Telescope image of NGC 7023, with the 
Spitzer/IRS SH field of view outlined by the black box. This image 
is composed from four filters: wide-band filters at 475 and 625 nm, 
a long-pass filter at 850 nm, and a narrow-ba nd filter at 658 nm . 
The exciting star, HD 200775, a B2.5 Ve star (|Finkenzeller|1985|), 
is indicated by the black cross towards the lower lett. 

irsclean, provided by the SSG. Spectra were extracted 
using the SSG pipeline modules profile, ridge, and 
extract. Extracted spectra were coadded to produce 
one spectrum per nod position, per instrument module. 
The products were then calibrated and stitched together 
to produce the final spectra. 

The SINGS data we re also provided in reduced form 
by Smith et al. (2007b). These data were processed with 
IRS pipeline version Sl 4 and reduced with the CUBISM 
tool (Smith et al.||2007a), including scaling and stitching 
of spectral segments to construct the final low-resolution 
spectra. 

The SH maps of NGG 7023, Girrus lb and Girrus 3 
were reduced in a similar manner to that of NGG 2023, 
which was described in paper 1. The raw data were pro¬ 
cessed by the SSG, using pipeline version S18.7. The 


reduction pr ocess was completed using CUBISM (Smith 


et al. 


mg. 


2007a|), including co-addition and bad pixel clean- 
'he extraction for NGG 7023 was accomplished by 


stepping a 2 X 2-pixel aperture across the map. 

3. ANALYSIS 
3.1. Spectral inventory 

A typical IRS low-resolution spectrum is shown in Fig¬ 
ure represented by the Hll region MSX LMG 1306. A 
rising mid-to-far-IR dust continuum is common in many 
objects in the sample. There are also several atomic 
emission lines frequently present, including [Ar ll] 6.99 
pm, [Ar hi] 8.99 pm, [S iv] 10.51 pm, [Ne ll] 12.81 pm, 
[Ne hi] 15.56 pm, and [S Hi] 18.71 pm. We frequently 
observe the S(3) 9.7 pm, S(2) 12.3 pm and S(l) 17.0 
pm rotational lines of H 2 , and occasionally G^o fullerene 


pm rotational lines or JT 2 , ana occasionally Gan miierene 
emission at 17.4 and 18.9 pm (|Gami et al.|2010t |Sellgren 
et al.|2010|). The PAH features at 6.2, 7.7, 8.6, 11.2, 12.7, 
15.8, 16.4, 17.4, and 17.8 pm are visible in all sources. 
Note that the 17.4 pm band arises from both PAH and 
Geo emission (Sellgren et al. 2010). In addition to the 
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Table 1 

Sample 


Object 

Type 

Resolv.“ 

SSID'’ 

AOR key 

Extent^ ["] 

RA [J2000] 

Dec. [J2000] 

HD269211 


Integ. 

4250 

19151360 

9.7 

05 12 30.26 

-70 24 21.75 

HD32364 

Ho?"* 

Integ. 

4122 

19150592 

11.2 

04 57 14.28 

-68 26 30.55 

IRAS05192-6824 

UGHiU 

Integ. 

4313 

19011584 

6.1 

05 19 06.84 

-68 21 36.41 

IRSX104 

Hii 

Integ. 

104 

22425088 

- 

05 24 13.30 

-68 29 58.98 

IRSX4087 

Hii 

Integ. 

4087 

11239168 

5.3 

04 55 06.47 

-69 17 08.30 

IRSX4713 

YSOf 

Integ. 

4713 

16943616 

6.3 

05 40 12.02 

-70 10 05.83 

IRSX4729 

YSQf 

Integ. 

4729 

16943616 

5.9 

05 40 46.83 

-70 11 22.73 

IRSX4755 

YSQf 

Integ. 

4755 

11239424 

7.0 

05 43 19.62 

-69 26 27.68 

MSXLMG1306 

Hii 

Integ. 

4059 

10958592 

6.5 

04 52 58.80 

-68 02 56.77 

MSXLMG217 

Hii 

Integ. 

4263 

10962176 

5.9 

05 13 24.68 

-69 10 48.12 

MSXLMG559 

Hii 

Integ. 

4419 

10963712 

- 

05 25 49.23 

-66 15 08.50 

MSXLMG764 

Hii 

Integ. 

4538 

10967040 

7.6 

05 32 52.62 

-69 46 22.84 

MSXLMG836 

Hii 

Integ. 

4527 

10966528 

6.4 

05 32 31.95 

-66 27 15.15 

MSXLMG889 

Hii 

Integ. 

4621 

10968576 

6.7 

05 38 31.64 

-69 02 14.92 

MSXLMG934 

Hii 

Integ. 

4652 

10969088 

13.4 

05 39 15.83 

-69 30 38.39 

N159-P2 

YSO 

Integ. 

4665 

12548352 

10.8 

05 39 41.88 

-69 46 11.94 

NGG 7023 

RN 

Extend. 

- 

3871232 

- 

21 01 32.50 

+68 10 18.0 

Girruslb 

ISM 

Extend. 

- 

4120832 

- 

17 32 53.80 

-33 10 01.4 

Girrus3 

ISM 

Extend. 

- 

4119296 

- 

16 03 38.70 

-52 18 01.7 


Notes. Additional sources 
“Whether the source 


(not listed) are from|Smith et al.| (|2007b|); |Boersma et al.l (|2010|); |Peeters et al.| 

, , :e is resolved or not (i.e., integrated spectrum or extende d ); mleigrateci source s are 

trom Surveying the Agents of Galaxy Evolution Spectroscopy program (SAGE-Spec; [Kemper et al.|2010| >, data 
release 3. ^SAGE-Spec ID number. “ Point sources are tho se with no extent listed. ^Object type uncertain in 
literature. ®Ultra-compact Hii region ( Beasley et al.|l996| . -^Glassification from Woods et al. (in prep.). 



Figure 2. A typical IRS low-resolution spectrum from the SAGE-Spec sample is shown, that of LMG Hll region MSX 1306. Prominent 
PAH features and atomic/molecular emission lines are visible (see labels at top). A local continuum is identified with a spline fit (red line). 
The continuum in the 5-10 pm range can be fit with or without a continuum point between the 7.7 and 8.6 pm features (see dashed red 
line). Three plateaus are identified underneath the PAH emission bands (green shaded region). 
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aforementioned PAH emission bands, there are several 
weaker PAH bands visible at 5.25, 5.75, 6.0, 11.0, 12.0, 
13.5 and 14.2 jam. Strong emission plateaus from 5-10 
pm, 10-15 pm, and 15-18 pm are observed. Additionally, 
som e spectra display a strong silicon carbide feature at 1 1 
pm ( Treffers fc Cohen||1974 Bernard-Salas et ar]|2009 ). 


3.2. Spectral decomposition 

A local spline continuum is fit to each spectrum to 
isolate the PAH emission features (Figure [^, a m ethod 
that has been used in several previous studies (e.g., Hon 


et al. 12001] IPeeters et al.|20021|Bernard-Salas et al. ^ 


)ny 


et al. |2010t|Stock et al.|2013||2014||. Spl: 

its are chosen to be aaiacent to emissi( 


Boersma ( 


line an- 


chor points are chosen to be adjacent to emission fea¬ 
tures, and are allowed to move within a predefined win¬ 
dow. When necessary, smoothing and/or a process that 
locates local minima are used to improve anchor posi¬ 
tions. In addition to the local continuum, a global con¬ 
tinuum is defined. This continuum underscores the emis¬ 
sion plateaus at 5-10 pm, 10-15 pm, and 15-18 pm. The 

7.7 and 8.6 pm features can be drawn with or without 
a spline anchor point between them, depending on the 
decomposition. Here, we use such an anchor point near 
8.2 pm. The molecular hydrogen and atomic emission 
lines wi ll be dealt with separately and are discussed in 
Section 13.31 

Other methods for measuring the PAH bands are pos¬ 
sible and are considered in the literature. These include 


PAHEIT tool; Smith et al. 

2007b) or Lorentzian 

files (Boulanger et al.| 1998 

Galliano et al. 

2008b 

has been shown that the decomposition metl 

rod cn 


It 


affects the measured band intensities, but it does not 
affect the t rends found in the correlations of band inten- 
sity ratios ( Smith et al.|2QQ7b Galliano et al.j2QQ8 ). The 
spline decomposition niethod is chosen here lor compar¬ 
ison with previous 1 5-20 pm PAH studies (paper I, II, 
Boersma et al.|[2MQ ). 


3.3. Measurement methods 

The plateaus are measured by direct integration be¬ 
tween the global continuum and local continuum (see 
Figure]^. After subtracting the local continuum, PAH 
band fluxes in the 5-15 pm range are measured by direct 
integration, except for the 6.0, 11.0 and 12.7 pm bands. 
The 6.0 and 11.0 pm bands are blended with the 6.2 
and 11.2 pm bands in low-resolution spectra, respec¬ 
tively. Therefore, Gaussians are used to separate these 
subfeatures and determine their fluxes (see paper II for 
details). Special care is required when measuring the 

12.7 pm PAH feature, as it blends with the 12.8 pm [Ne 
ll] line and the 12.3 pm H 2 line (when present). We 
separate these components by scaling the av erage 12.7 
pm PAH band profile (from |Hony et al.||2001l) to fit the 
data between 12.45 — 12.6 pm, where only PAH emission 
is expected (Figure]^. After subtracting the scaled PAH 
profile, the remaining lines are fit with Gaussians. This 
method is identic al to the decomposition described by 
Stock et al.|(|2014|), though we use the range 12.45 — 12.6 
pm instead of 12.4 — 12.6 pm. This is because the 12.3 
pm H 2 line is strong enough in some of our spectra to 
produce non-negligible (albeit minor) emission redward 
of 12.4 pm. We test for systematic effects in this de¬ 
composition by constructing an artificial spectrum. It is 


composed of the av erage 12.7 pm PAH band profile from 


Hony et al. (2001|) plus two Gaussians, for the 12.3 H 2 
line and the 12.8 [Ne II] line, and a variable amount of 
noise. We find that the derived PAH flux corresponds 
to the input artificial PAH flux within 2 to 3% when the 
integrated strength of the [Ne II] line is less than ten 
times that of the 12.7 PAH band. This is true for all 
our sources. The 12.3 pm H 2 line in our sample does not 
influence the derived PAH flux when decomposed this 
way. 

The 15-20 pm PAH bands and the atomic/molecular 
lines are also measured through Gaussian decomposition. 
The 16.4 pm band is slightly asymmetric, while the other 
15-20 pm PAH bands are symmetric (paper I). There is 
frequent blending between the 15.55 pm [Ne ill] line and 
the 15.8 pm PAH feature, as well as between the 17.1 
pm H 2 line and the 17.4 pm PAH band in low-resolution 
spectra. Emission by the fullerene Geo is k nown to be 


present at 17.4 and 18.9 pm in some objects (Garni et al 
mollSellgren et al.||2010| ). The 17.4 ^In Geo emission is 


approximately half of the intensity of the 18.9 pm band 


(estimated to be ^0.4 in Sellgren et al. 


2010, ^0.6 in 


~] 2010| and ^0.5 in |J3ernard-balas et n?]|2012 
1.5). The 18.9 pm band is observed in some o: 


Garni et al 
we adopt 0.5). The 18.9 pm band is observed in some of 
our spectra, typically blended to some degree with the 

18.7 [S III] line in low-resolution spectra. We examine 
the presence of Geo emission by fi tting the data with the 
model of Geo emission from TCI ( Cami et al.||2010 ) and 
an appropriate Gaussian for the [S ill] 18.71 pm line. In 
the sample of integrated sources, we detect Geo emission 
in four objects: IRSX 4755 (with associated signal-to- 
noise ratio, or SNR, of 3.4), IRSX 4087 (5.6), HD 32364 
(5.5) and HD 269211 (5.9). We correct our measured 17.4 
pm band flux by subtracting the expected 17.4 pm Geo 
component when the 18.9 pm Coo band is present. The 
residual flux, if significant, is then attributed to the 17.4 
pm PAH band. 

The measured PAH and Coo fluxes of the SAGE-Spec 
sources and the two cirrus sources are presented in Ta- 
blel The fluxes of the SINGS sources are shown in 
Tabled 


4. RESULTS 

We examine correlations between the 15-20 pm bands 
and the spatial distribution of the PAH emission in the 
extended RNe. NGC 2023S and NGC 2023N refer to the 
south and north maps of NGC 2023, respectively (papers 
I and II). Band fluxes in the correlation study are typi¬ 
cally normalized to the 11.2 pm band. Normalization is 
required because PAH abundances and object distances 
will vary from source to source. Only band ratios with 
3cr detections are included. The weighted Pearson cor¬ 
relation coefficients for all correlations are presented in 
Tabled 


4.1. Correlation plots 
4.1.1. The 164, M4 and 17.8 pm hands 

Correlation plots of the 16.4, 17.4, and 17.8 pm bands 
are presented in Eigured We address these in turn. 

17.8 vs. 16.4 V^m. These band fluxes are the strongest 
correlated for the integrated sources (correlation coef¬ 
ficient r = 0.93). Some segregation by object type is 
observed: two evolved stars sit at the lowest ratios of 
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Table 2 

PAH emission band fluxes: SAGE and ISM sources 


Object 

6.2 

|j.m 

7.7 

|j.m 

8.6 

|j.m 

11.2 

1 |j.m 

12.7 

|j.m 

15.8 |J.m 

16.4 ]xm 

17.4 [im 

17.8 ]xm 

Plateau 

18.9 |j.m Ceo 

HD269211 


- 


- 

3.64 

(1.14) 

13.09 

( 2 . 02 ) 

22.66 

( 1 . 51 ) 


- 

1.77 

( 0 . 11 ) 

0.42 

(0.11) 

0.56 

(0.10) 

8.81 

(0.88) 

1.07 

(0.18) 

HD32364 


- 


- 


- 

24.65 

( 1 . 68 ) 

23.01 

( 1 . 04 ) 


- 

2.02 

( 0 . 12 ) 

0.42 

(0.08) 

0.79 

(0.14) 

16.91 

(1.69) 

1.43 

(0.26) 

IRAS05192-6824 

14.10 

(0.77) 

28.35 

( 1 . 85 ) 

4.44 

(0.98) 

14.75 

( 1 . 70 ) 

14.35 

( 1 . 56 ) 


- 

1.54 

( 0 . 19 ) 


- 



7.99 

(0.80) 


- 

IRSX104 

8.33 

(1.53) 

10.40 

( 1 . 15 ) 


- 

6.45 

( 0 . 71 ) 

3.28 

( 0 . 75 ) 

0.44 

(0.12) 

0.72 

( 0 . 14 ) 


- 


- 


- 


- 

IRSX4087 

9.41 

(1.15) 

14.48 

( 1 . 34 ) 

3.50 

(0.75) 

10.60 

( 0 . 67 ) 

10.07 

( 0 . 74 ) 

0.63 

(0.13) 

1.13 

( 0 . 15 ) 


- 

0.48 

(0.12) 

5.27 

(0.53) 

1.62 

(0.29) 

IRSX4713 

5.76 

(0.64) 

8.47 

( 0 . 60 ) 

2.05 

(0.20) 

5.43 

( 0 . 18 ) 

1.48 

( 0 . 18 ) 



0.41 

( 0 . 10 ) 


- 


- 


- 


- 

IRSX4729 

6.95 

(1.00) 

10.65 

( 1 . 23 ) 

2.17 

(0.64) 

5.82 

( 0 . 57 ) 

1.92 

( 0 . 29 ) 


- 

0.41 

( 0 . 08 ) 


- 


- 

2.04 

(0.20) 


- 

IRSX4755 

18.30 

(2.66) 

23.74 

( 3 . 48 ) 

5.13 

(0.94) 

10.38 

( 1 . 17 ) 

6.02 

( 0 . 91 ) 


- 

0.82 

( 0 . 11 ) 


- 


- 


- 

0.34 

(0.10) 

MSXLMC1306 

21.60 

(5.09) 

38.72 

( 5 . 16 ) 

6.76 

(1.38) 

20.59 

( 1 . 33 ) 

12.46 

( 0 . 86 ) 

0.48 

(0.10) 

1.64 

( 0 . 10 ) 


- 

0.46 

(0.09) 

8.48 

(0.85) 


- 

MSXLMC217 

27.09 

(4.94) 

52.02 

( 4 . 33 ) 

7.27 

(2.00) 

27.14 

( 1 . 43 ) 

24.50 

( 1 . 25 ) 


- 

2.30 

( 0 . 25 ) 


- 


- 

11.70 

(1.17) 


- 

MSXLMC222 


- 


- 

3.10 

(0.33) 

6.97 

( 0 . 31 ) 


- 


- 

0.82 

( 0 . 11 ) 


- 


- 


- 


- 

MSXLMC559 

31.00 

(1.43) 

56.55 

( 1 . 39 ) 

7.92 

(1.01) 

18.46 

( 0 . 87 ) 

17.97 

( 1 . 30 ) 


- 

1.86 

( 0 . 61 ) 


- 


- 


- 


- 

MSXLMC764 


- 

76.81 (13.82) 


- 

39.66 

( 5 . 74 ) 

33.30 

( 3 . 61 ) 

1.85 

(0.33) 

4.40 

( 0 . 18 ) 

0.96 

(0.11) 

2.09 

(0.31) 


- 


- 

MSXLMC836 

22.14 

(0.67) 

34.98 

(1.29) 

5.80 

(0.52) 

19.88 

( 0 . 57 ) 

9.74 

( 0 . 41 ) 


- 

1.28 

( 0 . 12 ) 



0.52 

(0.13) 

7.88 

(0.79) 


- 

MSXLMC889 

61.10 

(7.74) 

154.46 

(6.41) 

20.51 

(4.29) 

70.84 

( 5 . 38 ) 

80.30 

( 3 . 13 ) 


- 

8.00 

( 0 . 17 ) 

1.47 

(0.28) 

2.79 

(0.34) 

36.50 

(3.65) 


- 

MSXLMC934 


- 


- 


- 

50.93 

( 1 . 71 ) 

52.89 

( 2 . 34 ) 


- 

4.88 

( 0 . 54 ) 

0.97 

(0.32) 

1.94 

(0.47) 

23.38 

(2.34) 


- 

N159-P2 


- 


- 


- 

32.70 

( 0 . 85 ) 

21.82 

( 0 . 99 ) 


- 

5.01 

( 0 . 80 ) 



1.93 

(0.63) 


- 


- 

Cirruslb 


- 


- 


- 

52.84 

( 1 . 08 ) 

33.84 

( 0 . 87 ) 

1.60 

(0.29) 

3.60 

( 0 . 48 ) 

1.37 

(0.28) 

1.30 

(0.32) 


- 


- 

Cirrus3 


- 


- 


- 

49.67 

( 0 . 94 ) 

27.65 

( 0 . 61 ) 

1.91 

(0.18) 

3.38 

( 0 . 18 ) 

1.22 

(0.23) 

0.95 

(0.17) 


- 


- 

Notes. In units of 10 
Uncertainties are given in 

W/m^, except the two extended cirrus i 
parentheses. 

sources, which 

are in 

units of 10 ^ 

W/m' 

^sr. 





Table 3 

PAH emission band fluxes: SINGS sources 

Object 

6.2 

pm 

7.7 pm 

8.6 

pm 

11.2 pm 

12.7 pm 

15.8 pm 

16.4 pm 

17.4 pm 

17.8 pm 

Plateau 

MRK33 

7.82 

(1.26) 

13.88 

(1.49) 

3.33 

(0.39) 

5.43 (0.30) 

3.72 (0.29) 

- 

0.45 

(0.09) 

0.16 

(0.04) 

0.16 (0.05) 


- 

NGC0337 


- 

17.62 

(1.62) 

3.58 

(0.45) 

7.84 (0.38) 

3.83 (0.30) 

- 

0.56 

(0.10) 


- 

- 

3.47 

(0.35) 

NGC0628 


- 




- 

3.48 (0.51) 

- 

- 

0.25 

(0.05) 


- 

- 


- 

NGG0855 


- 


- 


- 

3.50 (0.90) 

- 

- 

0.31 

(0.08) 


- 

- 


- 

NGG0925 


- 


- 


- 

4.72 (0.92) 

- 

- 

0.27 

(0.06) 


- 

- 


- 

NGG1097 

73.72 

(1.36) 

157.03 

(1.58) 

24.29 

(0.47) 

65.07 (0.42) 

38.90 (0.38) 

1.07 (0.22) 

3.58 

(0.42) 

1.56 

(0.45) 

- 

25.74 

(2.57) 

NGG1482 

100.70 

(1.51) 

230.20 

(2.04) 

31.10 

(0.76) 

74.11 (0.71) 

- 


4.36 

(1.29) 


- 

- 



NGG1512 




- 


- 

9.05 (1.48) 

- 

- 

0.47 

(0.08) 


- 

- 

3.42 

(0.34) 

NGG2403 

7.51 

(1.29) 

10.42 

(1.37) 


- 

5.13 (0.62) 

- 

- 

0.43 

(0.05) 


- 

- 


- 

NGG2798 






- 

47.55 (0.90) 

29.26 (0.84) 

0.83 (0.23) 

2.07 

(0.35) 

1.10 

(0.30) 

- 


- 

NGG2976 


- 


- 


- 

4.56 (0.64) 

- 

- 

0.38 

(0.05) 


- 

- 

1.83 

(0.18) 

NGG3049 


- 


- 


- 

5.41 (0.75) 

3.25 (0.68) 

- 

0.28 

(0.06) 


- 

- 



NGG3184 


- 


- 


- 

8.38 (1.08) 

5.47 (1.02) 

- 

0.31 

(0.04) 


- 

- 


- 

NGG3265 


- 


- 


- 

8.15 (0.75) 

3.89 (0.74) 

- 

0.46 

(0.11) 


- 

- 


- 

NGG3351 

16.01 

(1.94) 

53.97 

(2.09) 

5.30 

(0.82) 

17.63 (0.84) 

11.64 (0.85) 


1.04 

(0.15) 


- 

- 

12.92 

(1.29) 

NGG3521 

13.44 

(1.60) 

43.05 

(1.72) 

6.01 

(0.59) 

16.56 (0.52) 

- 

0.33 (0.07) 

1.21 

(0.10) 

0.36 

(0.08) 

0.33 (0.10) 

13.29 

(1.33) 

NGG3621 

15.28 

(1.68) 


- 


- 

12.12 (0.63) 

4.96 (0.67) 

- 

0.74 

(0.06) 

0.21 

(0.06) 

- 

4.64 

(0.46) 

NGG3627 

16.27 

(1.86) 

27.32 

(2.04) 

4.40 

(0.71) 

19.59 (0.86) 

- 

0.44 (0.08) 

0.86 

(0.12) 


- 

- 

10.66 

(1.07) 

NGG3773 




- 



4.12 (0.89) 

- 

- 




- 

- 



NGG3938 


- 


- 


- 

3.00 (0.20) 

1.48 (0.20) 

- 

0.18 

(0.05) 


- 

- 

2.18 

(0.22) 

NGG4254 


- 


- 


- 

23.23 (1.35) 

- 

- 

1.61 

(0.16) 


- 

- 

8.40 

(0.84) 

NGG4321 

33.00 

(2.12) 

69.04 

(2.58) 

10.03 

(1.01) 

27.13 (1.07) 

14.83 (1.10) 

- 

1.26 

(0.15) 

0.54 

(0.14) 

- 

11.56 

(1.16) 

NGG4536 

91.38 

(2.19) 

192.04 

(2.54) 

27.99 

(0.79) 

66.01 (0.70) 

- 

1.13 (0.36) 

2.86 

(0.55) 



- 

29.12 

(2.91) 

NGG4559 






- 

5.29 (0.47) 

2.84 (0.41) 

- 

0.36 

(0.06) 

0.12 

(0.04) 

- 



NGG4569 


- 


- 

5.89 

(1.16) 

24.68 (1.24) 

- 

0.44 (0.09) 

1.23 

(0.09) 


- 

- 

14.15 

(1.41) 

NGG4579 


- 


- 



4.76 (1.12) 

- 

- 

0.40 

(0.11) 


- 

- 



NGG4631 


- 


- 


- 

30.72 (0.55) 

- 

0.47 (0.13) 

2.09 

(0.20) 


- 

- 


- 

NGG4736 

34.32 

(1.62) 

40.69 

(2.13) 

9.73 

(0.99) 

48.69 (1.08) 

- 

0.89 (0.12) 

2.03 

(0.18) 

0.47 

(0.14) 

0.51 (0.12) 

16.99 

(1.70) 

NGG4826 

44.40 

(2.21) 

108.11 

(2.58) 

18.10 

(1.32) 

42.38 (1.44) 

- 

0.93 (0.11) 

2.30 

(0.18) 

0.72 

(0.16) 

0.72 (0.10) 

19.41 

(1.94) 

NGG5033 

19.75 

(1.76) 

36.88 

(1.75) 

6.15 

(0.79) 

17.27 (0.81) 

8.47 (0.74) 

- 

0.91 

(0.14) 

0.33 

(0.10) 

0.29 (0.09) 

9.01 

(0.90) 

NGG5055 

16.65 

(1.42) 

44.23 

(1.76) 


- 

21.37 (0.81) 

- 

- 

0.97 

(0.10) 

0.30 

(0.08) 

0.30 (0.06) 

10.06 

(1.01) 

NGG5194 

27.32 

(1.27) 

56.63 

(1.59) 

10.35 

(0.89) 

29.35 (0.88) 

- 

0.48 (0.10) 

1.28 

(0.14) 

0.43 

(0.11) 

0.34 (0.09) 

12.83 

(1.28) 

NGG5195 



32.48 

(2.59) 


- 

36.03 (1.28) 

- 

- 

1.93 

(0.46) 



- 

15.19 

(1.52) 

NGG5713 

36.36 

(1.25) 

75.25 

(1.62) 

11.23 

(0.72) 

29.04 (0.69) 

- 

- 

1.78 

(0.16) 


- 

- 

12.49 

(1.25) 

NGG5866 




- 


- 

8.24 (0.86) 

- 

0.22 (0.05) 

0.34 

(0.07) 


- 

- 


- 

NGG6946 

90.92 

(1.93) 

189.33 

(2.10) 

29.22 

(1.11) 

61.51 (1.05) 

- 

- 

4.34 

(0.41) 


- 

- 

33.86 

(3.39) 

NGG7331 




- 



15.56 (1.26) 

- 

0.27 (0.08) 

0.83 

(0.08) 


- 

- 

7.93 

(0.79) 

NGG7552 

81.57 

(1.84) 

173.99 

(2.22) 

26.90 

(0.97) 

71.78 (0.96) 

47.74 (0.95) 

- 

4.79 

(0.97) 


- 

- 


- 


Notes. In units of 10 W/m^. Uncertainties are given in parentheses. No 3(7 detections of 18.9 pm Ceo emission were 

found for these data. 
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Wavelength [/im] 


Figure 3. An example decomposition of the 12.7 pm complex. 
The continuum-subtracted spectrum (solid line) is fit with the 
scaled m ean 12.7 urn PAH band profile (dark gray shading), as 
found by Hony et al. ]l |200l| l. The fit is performed between 12.45 
and 12.65 pm (denoted by the dotted vertical lines). The fitted 
PAH profile (dark gray shading) can then be subtracted from the 
spectrum, leaving a remainder of the 12.3 pm H 2 emission line and 
the 12.8 pm [Ne ll] line (light gray shading). Returning to the orig¬ 
inal continuum-subtracted data, subtracting the two atomic lines 
then leaves only the inferred PAH 12.7 pm emission (dashed line). 


Table 4 

Weighted Pearson correlation coefficients 


Correlation 

Integrated 

sources 

N7023“ 

N2023N* 

N2023S'’ 

17.4 vs. 16.4 

0.67 

0.80 

0.57 

0.45 

17.8 vs. 16.4 

0.93 

0.85 

0.72 

0.82 

17.8 vs. 17.4 

0.63 

0.53 

0.97 

0.47 

16.4 vs. 15.8 

0.74 

0.22 

- 

-0.28 

17.4 vs. 15.8 

0.58 

0.37 

- 

-0.14 

17.8 vs. 15.8 

0.75 

0.40 

- 

-0.08 

15.8 vs. 11.2 

-0.28 

0.38 

- 

0.89 

Plat. vs. 11.2 

0.42 

0.47 

0.92 

0.87 

Plat. vs. 15.8 

0.28 

0.46 

- 

0.88 

12.7 vs. 11.0 

0.51 

0.95 

0.88 

0.91 

16.4 vs. 11.0 

0.36 

0.92 

0.92 

0.95 

17.4 vs. 11.0 

-0.14 

0.84 

0.69 

0.51 

17.8 vs. 11.0 

-0.26 

0.87 

0.89 

0.81 

16.4 vs. 12.7 

0.28 [0.68]^ 

0.89 

0.88 

0.92 

17.4 vs. 12.7 

0.69 

0.83 

0.55 

0.35 

17.8 vs. 12.7 

0.23 [0.66]^ 

0.86 

0.53 

0.71 


Notes. “ See individual correlation plots for coeffi cients 
organized by region; regions defined in Figure [7o| ^ See also 
paper I, H. ^ Hll regions alone. 


16.4/11.2 and 17.8/11.2 (0.030 and 0.010, respectively), 
galaxies slightly higher at (0.060, 0.015), and Hll re¬ 
gions at the highest values (centered near 0.10, 0.035). 
The galaxies and H ll regions each span a range covering 
roughly a factor of two on each axis. The 17.8 and 16.4 
pm bands also correlate in the RNe. NGC 7023 has a 
correlation coefficient of r = 0.85 and spans nearly the 
same variation in 16.4/11.2 and 17.8/11.2 as the inte¬ 
grated sources. NGC 7023 has a cluster of data points 
centered near (0.13, 0.030), which we determine to orig¬ 
inate in the portion of the map closest to the exciting 
star (the upper third of the map; see Figure [Tq|). NGC 
2023S has a correlation coefficient of r = 0.82 and has a 


line of best fit close to that of NGC 7023, despite span¬ 
ning half the range in abscissa; NGC 2023N has a similar 
range and correlation coefficient (r = 0.72). Its line of 
best fit is closer to the fit for the integrated sources, but 
it should be noted it also has far fewer pixels than those 
contained in the NGC 7023 or NGC 2023S maps given 
the SNR criterion. 

17.4 vs. 16.4 [im. The 17.4 and 16.4 pm band fluxes 
are observed to correlate. NGC 7023 exhibits a correla¬ 
tion coefficient of r = 0.80, while NGC 2023S and 2023N 
have associated coefficients of r = 0.45 and r = 0.57, re¬ 
spectively. NGC 7023 and NGC 2023 have similar lines 
of best fit. The integrated sources however display a 
trend on a much steeper gradient, aside from two lone 
H ll regions, and have a correlation coefficient of r = 0.67. 
The integrated sources and NGC 2023 generally span the 
same abscissa and ordinate ranges, but NGC 7023 is set 
to higher ratios of 16.4/11.2, with only modest overlap 
with NGC 2023 and the integrat ed sources. This correla - 
tion was previously observed by Berne & Tielens (2012). 

17.8 vs. 17.4 [im. The 17.8 and 17.4 pm bands also 
correlate. The integrated sources and NGC 7023 span 
similar ranges in abscissa and ordinate values, with as¬ 
sociated coefficients r = 0.63 and r = 0.53, respectively. 
NGC 2023 spans half the range of 17.8/11.2 in compar¬ 
ison, though a similar extent is observed in 17.4/11.2. 
NGC 2023S has correlation coefficient of r = 0.47 and 
NGC 2023N of r = 0.97, though the latter has far fewer 
data points included (12 versus 52, respectively). The 
individual sources again fall upon a slightly different 
gradient, as they appear to be displaced toward higher 
17.4/11.2 values. The outlying Hll regions in this figure 
are also the outliers in the correlation plot of the 17.4 
and 16.4 pm bands. Systematic problems may affect the 
17.4 pm band measurement in our sample as its derived 
fiux is dependent on an accurate measurement of the 18.9 
pm Geo band. When the latter cannot be observed (or 
reliably measured) it is possible that the derived 17.4 
pm PAH band flux is overestimated, as the 17.4 pm Geo 
contribution cannot be removed, if present. This issue is 
particularly prominent in low-resolution spectra, which 
comprise most of the integrated sources in this sample 
(apart from the two ISM cirrus sources), as the 18.9 
pm Geo band blends with the [S Hi] 18.71 pm line. The 
RNe were observed in high resolution and are less suscep¬ 
tible to this type of overestimation, which may explain 
why the gradients observed in the correlation plots in¬ 
volving the 17.4 pm band differ for integrated sources 
versus RNe (e.g.. Figure®. 

Correlations with the lz.7 pm band. The 16.4, 17.4 and 
17.8 pm bands are observed to correlate with the 12.7 
pm band (Figure [^. NGC 7023 displays a high correla¬ 
tion coefficient (r = 0.89), though there is deviation from 
the linear fit in our data at high abscissa values. In our 
data, there is a cluster of data points at 12.7/11.2=0.75, 
16.4/11.2=0.13, all of which originate in the third of the 
map nearest the exciting star in NGC 7023. Regardless, 
we find a strong correlation between these band ratios 
in NGC 7023, peaking in the middle portion of the map 
(r = 0.94; see Figure [1Q| for map region definitions). The 
gradient of the NGC ^023 data is consistent with those 
of NGC 2023 (both north and south). The integrated 
sources of this study show no trend as a whole, but they 
do tend to cluster about a position coincident with the 
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16.4/11.2 


c\j 


00 


(b) • NGC 7023 near (r = 0.63) 



16.4/11.2 



Figure 4. Correlation plots of flux ratios amongst the 16.4, 17.4 
and 17.8 |xm features. The fllled blue and orange density con¬ 
tours represent the number densities of measurements of NGC 2023 
North and South, respectively (papers I and II). Each displays ten 
linear contour levels, computed with a Gaussian kernel. The data 
points with green errorbars are fro m NGC 7023, separated into 
three regions (deflned in Figure [To| : data nearest to the exciting 
star (red dots); data from the center of the map, the peak PAH 
emission region (black dots); and data furthest from the exciting 
star (white dots). Each grouping is shown with an associated line 
of best fit when the correlation coefficient exceeds r = 0.50: inte¬ 
grated sources (black line), NGC 7023 (green line), NGC 2023N 
(blue line), NGC 2023S (orange line). 


center of the NGC 2023 data. The subset of Hll regions 
alone however do display a correlation (r = 0.68). Since 
the 16.4, 17.4 and 17.8 pm bands correlate (Figure®, 
it is not surprising that we find that the 17.4 and 1778 
pm bands correlate with the 12.7 pm band in RNe. We 
additionally observe that Hll regions display a correla¬ 
tion between the 17.8 and 12.7 pm bands. Nothing can 
be determined about Hll regions amongst the 17.4 and 
12.7 pm bands as there are too few sources meeting the 
detection criteria. The 12.7 pm band was previously ob¬ 
served to correlate with the 16.4 pm band in NGC 2023 
(paper I; r = 0.92 in the souther n map, r = 0.88 in the 
northern map) and in NGC 7023 (Boersma et al.||2014a 
their Figure 9). 

Correlation with the 6.2 pm band. We examine the 
flux ratios of the 16.4 and 6.2 pm bands in Figure® for 
the integrated sources in our sample. The 6.2 pm band 
cannot be measured in the SH RNe data as the wave¬ 
length cover age begins at ^ 10 pm. Data from [Boersma 


et al. (2010) are included, where we have replaced the 


authors' average SINGS galaxy data point with our indi¬ 
vidual SINGS sources. We observe no linear correlation 
(r= 0.11), though we note the lower-right quadrant is 
devoid of sources. Furthermore, we identify the Hll re¬ 
gions as occupying the upper-left quadrant of this fig¬ 
ure. If these are excluded from the data set, a mod¬ 
est correlation is observed (r= 0.47). A mostly linear 
trend (r^ = 0.76) bet ween these bands has b een identi¬ 
fied in NGC 7023 by|Boersma et al.| (|2014a|), their Fig¬ 
ure 10, albeit with some residual curvature. Our mea¬ 
surements are consistent with their measurements at low 
abscissa/ordinate values. 

Since we detect a correlation between the 16.4 and 12.7 
pm bands, and essentially none between the 16.4 and 
6.2 pm bands, we should expect no correlation between 
the 6.2 and 12.7 pm bands for consistency. We indeed 
detect no such correlation in our data alone. However, 
the 6.2 and 12.7 pm PAH bands have been observed to 
correlate for a sample of inte grated sources of various 
objec t types (iHony et al. 200 1) and within the RNe NGC 
7023 (jBoersnia et al.||2ol4ap and NGC 2023 (paper II). 
The lack of such a correlation in our sample of integrated 
sources is likely due to the small range of ratios probed 
by our data, as they span a small range of ratios relative 


to 


Hony et al. (2001). 

Correlations with the 11.0 pm band. We observe that 
the 16.4, 17.4 and 17.8 pm bands correlate with the 11.0 
pm band in NGC 2023 and NGC 7023 (Figure®. The 
NGC 7023 pixels furthest from and nearest to me star 
tend to individually cluster, while the region containing 
the peak of the PAH emission shows linear correlations 
between band intensity ratios (regions defined in Fig¬ 


ure 10). Th is clustering has been pr eviously observed, for 
example in Boersma et al. (2014a). The 12.7 pm band 
is also included for comparison, and the RNe display 
strong correlations between the 12.7 and 11.0 pm bands 
(r > 0.90). The integrated sources show no clear cor¬ 
relations between the 11.0 pm band and the 12.7, 16.4, 
17.4 and 17.8 pm bands, as they tend to cluster near 
11.0/11.2^0.060. They tend to be consistent with the 
observed RNe correlations for the 12.7 and 16.4 versus 
11.0 pm bands, but they deviate for the 17.8 pm and 
especially the 17.4 pm band. 
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12.7/11.2 


0.10 


0.08 

C\1 

0.06 

00 

1 ^ 

0.04 


0.02 


0.00 

0 


Figure 5. Correlation plots of the 16.4, 17.4 and 17.8 |xm band 
fluxes versus the 12.7 |xm band flux, each normalized to the 11.2 
|xm band flux. 


NGC 7023 mid (r = 0.83) 

NGC 7023 far (r = 0.93) 

■ NGC 7023 (r = 0.86) 

NGC 2023N (14 px, r = 0.53) 
NGC 20238(103 px, r = 0.71) 
' Integrated sources (r = 0.23) 
Hll regions (r = 0.66) 



0.4 0.6 

12.7/11.2 


4.1.2. The 15.8 \im band 

The 16.4, 17.4 and 17.8 jam bands do not appear 
to correlate with the 15.8 pm band (Figure]^. NGC 
2023S shows no correlation between these bands (e.g., 
r = —0.28 when comparing the 15.8 and 16.4 pm bands) 
and NGC 2023N has too few detections to be represented 
here (papers I and II). NGC 7023 generally appears to 
have a similar distribution to NGC 2023S, with no over¬ 
all correlation (e.g., r = 0.22). However, NGC 7023 
does seem to exhibit a “C”-shaped distribution in the 
16.4 versus 15.8 pm band intensity correlation, and to 
a lesser degree in the 17.8 versus 15.8 pm correlation 
plot. This may be indicative of a relationship present at 
16.4/11.2 ratios greater than 0.12, and 17.8/11.2 ratios 
greater than 0.03, but it is not conclusive. The inte¬ 
grated sources are shown to o ccupy three regi o ns of the 
correlation plots: the data of Boersma et al. (2010) at 
low abscissa values; the SINGb galaxies at an interme¬ 
diate position; and H II regions or cirrus sources towards 
the upper-right quadrant. Together these produce a cor¬ 
relation (with coefficient r = 0.74) but this is possibly 
spurious. Specifically, in the previous plots (Figure®, 
the identified correlations can be seen both individu- 
al ly within our sample and individually in the sample 
of Boersma et al. (2010). In contrast, these correlations 
involving the 15.8 pm band are not present in either sam¬ 
ple in isolation. Additionally, the closeness of the [Ne ill] 
15.56 pm line to the 15.8 pm PAH band in low-resolution 
spectra affects how well the PAH band can be measured. 
This is compounded because the 15.8 pm PAH band is 
weak and is near the inflection point of the rising 15-18 
pm continuum; the error on the continuum is not taken 
into account in the 3cr cut-off. Finally, the Pearson cor¬ 
relation coeffivinn^j^ known to be sensitive to outlying 
points (jDevlin et al. 1975), which is why the measured 
coefficients here are not necessarily meaningful, unlike in 
Figure 

The 15.8 and 11.2 pm bands were shown to correlate in 
NGG 2023 (paper I). We revisit this relationship in Fig¬ 
ure |^a). The data near the peak PAH emission re gion 
(i.e., m the middle of the NGG 7023 map, see Figure 10) 
are shown to correlate slightly (r = 0.46) and exhibit 
a gradient consistent with that in NGC 2023S. No cor¬ 
relation is identified for the integrated sources, though 
Hll regions tend to lie at relatively high 15.8/12.7 ra¬ 
tios when compared to evolved stars or galaxies. The 
derived fluxes of the 15.8 and 12.7 pm bands are in¬ 
fluenced by the presence of the [Ne II] 12.81 pm and 
[Ne III] 15.56 pm emission lines. These lines are present 
in Hll regions and PNe and may influence or explain 
the observed discrepancy. Figures ®b) and (c) com¬ 
pare the emission of the 15-18 pm gateau to those of 
the 11.2 and 15.8 pm bands, respectively, which were 
shown to correlate in NGC 2023 (paper I). We find that 
NGC 7023 has similar correlations within its peak emis¬ 
sion region; the other NGC 7023 regions do not show a 
clear correlation by themselves but they are consistent 
with the correlation seen in the NGC 7023 peak (“mid”) 
and NGC 2023. The integrated sources show no linear 
trend but there is significant clustering near the position 
(Plat./12.7, 11.2/12.7) - (1.5,3.5). 


4.2. Spatial distribution of 15-20 pm PAH emission 
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Figure 6. Correlation plots of the 12.7, 16.4, 17.4 and 17.8 pm bands versus the 11.0 pm band. 

4 . 2 . 1 . Map morphology 

We have thus far focused solely on correlations between 
bands. Paper II finds that the tightness of a correlation 
between any two bands is related to their spatial coin¬ 
cidence. However, even two bands with different spa¬ 
tial distributions could retain high correlation coefficients 
(e.g., the 8.6 and 7.7 pm bands, with r = 0.936; paper H). 
It is therefore important to determine if the correlations 
we identify are corroborated by their spatial morpholo¬ 
gies. To do this, we examine the spectral map of NGC 
7023 (Figure [1Q|) and compare it with the spectral map of 
NGC 2023 (papers I and II). The 16.4 pm data from this 
spectral cub e has been previously examined by [Sellgr"^ 



et al. 


(|2007|, and th e 11.0, 11.2, 12.7 , 16.4 pm bands and 
plateau emission by Boersma et al.|(2013, 2014a, 2015). 
Additionally, zones parallel to theP DK f ront have been 


6 . 2 / 11.2 

Figure 7. No linear relationships is observed between the 16.4 
and 6.2 pm bands (r = 0.11); however, the lower-right quadrant is 
conspicuously devoid of sources, and our data span a small range 
of values (c.f. [Boersma et al.|2014a| . 


identifi ed by Berne & Tielens (2012) and Boersma et al. 


(2014a); we introduce three such zones in Figure" 10 


The spatial distributions of the 12.7, 16.4 ana 17.8 
pm bands are very similar in appearance. Each of these 
is seen to be displaced slightly toward the irradiating 
star relative to the 11.2 pm band. The 17.8 pm band is 
slightly more compact than the 12.7 and 16.4 pm bands; 
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15.8/11.2 




Figure 8. Correlation plots between the 11.2, 15.8 |xm bands 
and the 15-20 |xm plateau (panels a, b, c, respectively). 





15.8/12.7 


Figure 9. Correlation plots involving the 11.2 and 15.8 pm bands 
and the 15-18 pm plateau, (a) The 15.8 and 11.2 pm bands, nor¬ 
malized to the 12.7 pm PAH flux, correlate in RNe. (b), (c) 15-20 
pm plateau correlations with the 11.2 and 15.8 pm bands in RNe. 
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Figure 10. Maps of the normalized PAH band fluxes towards NGC 7023 (13 x 15 pixels). Maps have been rotated 20° such that the star 
is at the top of each panel in this orientation. The original held of view is as indicated in Figure^ The black contours are constructed 
from the 11.2 pm band flux, normalized to unity. Bands are organized into vertical groupings (columns) based on similar characteri stic s 
(see text). The vertical dashed lines (shown only on the 11.2 pm map) denote the region used for the 5 pixels-wide radial cut (Figure p^ . 
We introduce in the lower-right panel three simple zones to distinguish regions of NGC 7023 data present in the correlation plots. 


we note it also has a limited SNR in comparison. In 
NGC 2023, the 12.7 and 16.4 pm bands are also nearly 
identical in appearance, and the 17.8 pm band is mildly 
more compact. The 17.4 pm band has some apparent 
differences from the 12.7, 16.4 and 17.8 pm bands in the 
spatial maps, despite correlating with these bands: its 
emission is more compact and appears to peak 10" closer 
to the exciting star in both NGC 7023 and NGC 2023, 
though the 17.4 pm band has a limited SNR in NGC 
7023. The 11.0 pm band is also displaced closer to the 
exciting star to a similar degree in both RNe. The 17.4 
pm band appears identical to the 11.0 pm band in NGC 
2023 prior to removing the Geo component; after correc¬ 
tion, the 17.4 pm band is slightly more compact than the 
11.0 pm band. 

The 11.2 pm band and the 15-18 pm plateau have 


coincident peak positions in both RNe, while the 15.8 
pm band peaks slightly further distant from the star. 
The 15.8 pm band is much less extended than the 11.2 
pm band and the plateau in NGC 7023, but it is only re¬ 
liably measured in a small portion of the map due to its 
low SNR. In NGC 2023, the 15.8 pm band has a higher 
SNR and it shows a spatial distribution extremely close 
to that of the 11.2 pm band. The 15-18 pm plateau in 
NGC 2023 is much broader in extent. 

4.2.2. PAH emission variability with projected distance 

In order to investigate the change in PAH emission 
with projected distance from the illuminating star, we 
introduce radial slices or cuts across the maps of NGC 
7023 and NGC 2023 (Figures [To] and [IT]) and analyze the 
individual band intensities as a function of the strength 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 


Figure 11. A normalized 11.2 |xm PAH intensity map of NGC 
2023 S (papers I and II). These data have been rotated 35 degrees 
counter-clockwise to be roughly perpendicular to the irradiating 
star’s radial vector. The outlined regions identify the 5 pixels- 
wide apertures used to create the radial slices. The exciting star is 
located upward in this projection, with estimated distance to each 
aperture as shown. 


of the radiation field (Figures and [T^. 

NGC 7023 is ir radiated by HD 2() Q77Fi a star of spec¬ 
tral type B2.5 Ve (Fink enzeller|198^, of effe ctive temper¬ 
ature 17QQQi^±lQ00A ( [Baschek et al.|l982| ). We find the 
strength of the FUV radiation field, Go, to be 1.2 x 10^ at 
12" and 1.4 x 10^ at 35" from HD 200775 F] This is con¬ 
sistent with|Berne & Tielens| (|2012|), who adopted values 
of 1 ± 0.7 X 10"^ at 12" and 1 ± 0.7 x 10^ at 35". 

The slices of NGC 7023 show that the 11.2 and 15.8 
pm bands and the 15-18 pm plateau peak furthest from 
the star, at a projected distance, d, of 49", or Go ^ 
7.0 X 10^. The 12.7, 16.4 and 17.8 pm bands peak closer 
to the star, at d^ 39"—44", or Go ^ 11 x 10^ — 8.5 x 10^. 
The 17.4 pm band peaks at d^ 40"(Go ^ 1.2 x 10^) and 
the 11.0 pm band has a wide peak from d^ 30 — 39", 
corresponding to Go ^ 1.9 x 10^ — 1.1 x 10^. The 
17.4 pm band profile is much narrower than the 11.0 
pm emission, but it is aligned with the far edge of the 
11.0 pm band. Slices of the 16. 4 and 18.9 pm band s 
in NGC 7023 were presented by Sellgren et al. (2007). 
These authors observed the 16.4 pm emission to peak at 
approximately 36" distant from the star. We observe this 
peak to lie more distant, near 44", but we use a differ¬ 
ent decomposition method. Overall, the spatial profiles 
of the 12.7, 16.4 and 17.8 pm bands are broader and 
more symmetric than those of the other bands. We find 
that we can reproduce the profiles of the 12.7 and 16.4 
pm bands extremely well by adding the 11.0 and 11.2 
pm profiles and renormalizing the result (see Figure 12). 

We perform a similar analysis on NGC 2023S, wnas 
data have higher SNR and overall intensities than NGC 
2023N’s (Figure 13). NGC 2023 is illuminated by a 


® Go is the FUV flux between 6-13.6 eV, expressed in terms 
of the Habing held, 1.2 x 10~^ erg/s/ cm^/sr ( [Habing |l968| . We 


Tielens 


(2005[: Go = 625^^-, where L 


use the formulation of 
is the luminosity of the star, y is the luminosity fraction between 
6-13.6 eV and d is the distance from the star, x was determined 
by evaluating a blackbody of effective temperature corresponding 
to the appropriate stellar spectral type. 


Go 

50 X 10^ 25 X 10^ 15 X 10^ 10 x 10^ 5 x 10^ 



Figure 12. Normalized PAH intensity along a projected cut 
across NGG 7023 directed towards the irradiating source, extracted 
at 12 p ositions, using a 5 pixels wide extraction aperture (see Fig¬ 
ure II^. The colors indicate groupings of bands with similar peak 
positions and/or spatial profiles along the cut. 


B1.5V star ( Mookerjea et al.|2QQ9 ). We derive Go values 
in a similar manner as for JNGC 7023. The PAH transi¬ 
tions occur from Go ^ 1.5x10^—3.0x10^, in rough agree¬ 
ment with the estima tes of 10^ — 10^ in the literature, 
albeit on the low end (Sheffer et al. 2011 paper I). The 
spatial stratification seen in NGC 7023 is also present in 
NGC 2023S. However, the transition occurs at quite dif¬ 
ferent values of Go than NGC 7023. This suggests that 
our simple approximation of Go is not accurate enough 
and other parameters such as density, gas temperature 
and/or projection effects may play a role. Specifically, 
the 11.0 and 17.4 pm bands peak nearest the star, with 
the 11.0 pm band showing a broad peak from approxi¬ 
mately 80-95" distant (Go ^ 2.9 x 10^ — 2.1 x 10^), and 
the 17.4 pm band being more compact near a distance of 
84"(Go ^ 2.6 X 10^). The 11.2, 15.8 pm bands and the 
15-18 pm plateau all peak near d^lll"(Go ^ 1.5 x 10^), 
but we note the plateau is slightly broader. The 12.7 
and 16.4 pm bands are fiat from d^80—105" (Go ^ 
2.9x10^ — 1.7x10^), again broader than the other emis¬ 
sion bands, as in NGC 7023. The 17.8 pm band is similar 
to the 12.7 and 16.4 pm bands but it peaks closer to the 
star (80-95", Go ^ 2.9 x 10^ — 2.0 x 10^), with a pro¬ 
file near that of the 11.0 pm band. Again, the 12.7 and 
16.4 pm band profiles are well realized by summing and 
renormalizing the 11.0 and 11.2 pm bands, albeit with a 
slight overestimate near 110". 


5. DISGUSSION 
5.1. Interpretation 

We first summarize the observational results, and sub¬ 
sequently, we discuss the implications in terms of charge 
and molecular structure, based on comparison with theo¬ 
retical studies. The observational results consist of PAH 
band intensity correlations for integrated and extended 
sources, as well as spatial measurements of PAH emission 
in spectral maps of NGC 7023 and NGC 2023. 

Correlations. The 16.4, 17.4 and 17.8 pm band 
fluxes are inter-correlated in all environments. Further¬ 
more, each of these three bands correlates with the 12.7 
pm band in the RNe. Amongst the sample of integrated 
sources, the 16.4 and 17.8 pm bands correlate with the 
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Figure 13. Normalized PAH intensity along a projected cut 
across NGC 2023S directed towards the irradiating source, ex¬ 
tracted at 8 p ositions, using a 5 pixels wide extraction aperture 
(see Figure HH. The colors indicate groupings of bands with simi¬ 
lar peak positions and/or spatial profiles along the cut. 

12.7 jLLm band only in Hll regions. Together, these four 
bands (12.7, 16.4, 17.4 and 17.8 jam) correlate strongly 
with the 11.0 JLLm band in RNe. The 15.8 JLLm band 
does not correlate with any of these bands, leaving it 
as the odd one out. However, the 15.8 jnm band is inter- 
correlated in RNe with the 15-18 jam plateau and the 
11.2 pm band. Finally, the 16.4 and 6.2 pm bands show 
no overall linear correlation for the integrated sources, 
though if we exclude H ll regions from this analysis they 
show a modes t correlation, in line with r esults found in 
the literature (Boersma et al.|2010 2014a). Our data are 
also isolated to relatively low ordinate and abscissa val- 
ne s, falling upon the lowe r end of the relationship found 
by Boersma et al. (2014a). 

Spatial maps. Tne RNe spatial maps show a common 
morphology between the 12.7, 16.4 and 17.8 pm bands. 
Similarity is also observed between the 11.2 and 15.8 
pm bands and the 15-18 pm plateau. Lastly, the 11.0 
and 17.4 pm bands appear to peak in similar locations. 
We invoked radial cuts to isolate the spatial origin of 
each band to examine these similarities. The 11.0 and 

17.4 pm bands peak closest to the star, in similar loca¬ 
tions, though the 11.0 pm band’s spatial profile is wider. 
The 11.2, 15.8 pm bands and 15-18 pm plateau peak 
furthest from the star, sharing very similar profiles. The 
12.7 and 16.4 band profiles peak in an intermediate loca¬ 
tion, with very similar shapes. The 17.8 pm band follows 
the 12.7 and 16.4 pm bands in NGC 7023 but its pro¬ 
file is more similar to that of the 11.0 pm band in NGC 
2023S. We observe that the profiles of the 12.7 and 16.4 
pm bands bands (and the 17.8 pm band in NGC 7023) 
can be reproduced by averaging the spatial profiles of the 
11.0 and 11.2 pm bands. 

5 . 1 . 1 . Charge 

We first examine the observational results in the con¬ 
text of ionization states. The 11.0 and 11.2 pm PAH 
bands have been firmly assigned to cationic and neu¬ 
tral PAHs, respe ctively (Hudgins & Allamandola 1999; 
Hony et al.||2001). The similarities of the 15.8 pni banc, 
to the 11.2 pm band in the maps, cuts, and the fact 
that they correlate suggest that the 15.8 pm band origi¬ 


nates in neutral PAHs. Likewise, the observational data 
suggests the 17.4 pm band is coincident with the 11.0 
pm emission, with which it correlates; we therefore asso¬ 
ciate the 17.4 pm band with cationic PAHs. The 12.7, 

16.4 and 17.8 pm bands are self-similar in the maps and 
cuts, and have an extent that envelops both the 11.0 and 
11.2 pm bands. In conjunction with the correlation plots, 
these data suggest that the 12.7, 16.4 and 17.8 pm bands 
arise from both cationic and neutral PAHs. 

Paper I suggested that the 17.8 pm band arises from 
both neutral and ionized PAHs, which we here confirm. 
Furthermore, paper I concluded that the 16.4 pm band 
is from species cospatial with the cationic portion of the 
emitting PAH populatio n. This result was stro ngly sup¬ 
ported by the results of Boersma et al. (2014a); the au¬ 
thors however noted that residual curvature m their cor¬ 
relations plots indicated that other population changes 
are at play. Here we identify two components of the 

16.4 pm band: one cospatial with cationic PAHs and 
one cospatial with neutral PAHs. The 12.7 pm band was 
al so shown to be co i nciden t with PAH cations in paper 
1. Rosenberg et al. (2011) used blind signal separation 
to identity two coniponents of the 12.7 p m band, which 
they as sociate with neutrals and cations. [Boersma et al.| 
(2014a) supported their assignment; this work also sup¬ 


ports a combined neutral and cationic origin for the 12.7 
pm band. Paper I suggested the 17.4 pm bands arises 
from doubly ionized PAHs and/or a subset of the cationic 
PAH population. We attribute the 17.4 pm bands to 
cations generally, and not necessarily a particular subset. 
The discrepancy may be due to how each paper traces 
cations: we use the 11.0 pm band as a tracer of cations 
here, whereas paper I used the 6.2 and 7.7 pm bands 
as tracers. Paper II finds that there are two popula¬ 
tions contributing to the 7.7 pm band (and thus the 6.2 
pm band, due to its similar morphology), and so the 7.7 
pm band ma y not be an ideal trac er of PAH cations. 
Additionally, Boersma et al.j ([2015 ) show that the 7.7 
pm PAH band strength may not be a reliable tracer for 
PAH cations. We attribute the 15.8 pm band to neutral 
PAHs, in agreement with paper 1. 

Theoretical studies of PAH vibrational modes have 
given us clu es about the orijS!;ins of the 15-20 pm emis - 


sion bands (Boersma et al. 
These often make use of the 


2010 Ricca et al. 2010). 


ASA Ames PAH IR^bic- 


trosco 

pic Database 

!] 

([Bauschlicher et al. 2010 

et al. 

r—1 

o 

CN 

)■ Spf 

3C1 

lically, the 16.4 pm band 


Boersma 


to arise Ironi a mixture of out-of-plane and in-plane vi- 
bration s, particularly due to elongation and compression 


modes (Ricca et al. 2010). The 17.4 pm band originates 
in a mixture of C-O-C in-plane and out-of-plane modes 
and the 17.8 pm band is attributed to a mixture of G-H 
and G-G-G out-of-plane bending. The origin of the 15.8 
pm band is unclear, though emission at 15.4 pm by G-G- 
C in-p lane bending modes was identified by [Ricca et al. 
(2010); this may or may not correspond to the astro¬ 


nomical 15.8 pm band. Using these analyses, we assert 
that the intrinsic intensities of the 15.8, 16.4, 17.4 or 17.8 
pm bands do not vary by more than approximately a fac¬ 
tor of 2 between charge states (anions, neutrals, cations). 

^ http://WWW.astrochem.org/pahdb/ 
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This estimate is based on Figure 13 of Ricca et al. (2010), 
which presents average spectra for all FAHs in the Ames 
database containing more than 50 carbon at oms, sep¬ 


arated by io nization state, and Figure 15 of Boersma 


et al.l ( 2010 ), showing the intrinsic intensities tor differ 


ent charge states of C 96 H 24 . In general then, the 15-20 
pm emission bands produce intensity variations that are 
much more subdued than the variations observed in the 
3-12 pm ra nge, which span approxim ately one order of 


. OP-I 

magnitude (|Allamandola et al.|1999). Emission by an¬ 
ions, neutrals and cations i s expected for all four 15-20 
pm bands in the database (Ricca et al.||2010). The 16.4 
pm band is expected to be tw ice as strong for cations 
than for the other charge states ( Bauschlicher et al .|2010 


Ricca et al. 2010). The 17.8 pm band appears to have 
similar intensities in each charge state, while the 17.4 
pm band should be st rongest in anions by a factor of 
two (Ricca et al. 2010). These authors computed spec¬ 
tra for eight neutral FAHs and noted similar intensities 
between the 17.4 and 17.8 pm band, possibly implying a 
correlation; such a correlation is observed i n this work. 


Lastl y, emission was observed at 15.4 pm by (Ricca et al. 


2010 ), with cations exceeding the intensities ol neutrals 


and anions by a factor of two, but it is unknown if this 
emission corresponds to the astronomical band at 15.8 
pm. 


Setting aside anions due to their expected very low 
abundances, we note that the prediction of similar emis¬ 
sion intensities of the 17.8 pm band between cations and 
neutral FAHs is consistent with our observational re¬ 
sults. This can be seen in the spatial maps and cuts, 
in which the 17.8 pm band’s spatial profile encompasses 
the 11.0 and 11.2 pm bands, centered between them. 
The 16.4 pm band, which is predicted to be twice as 
strong in cations than neutrals, is however somewhat at 
odds with the observational data: we observe it has the 
same spatial profile as the 17.8 pm band, with no prefer¬ 
ence or tendency towards the 11.0 pm emission over the 
11.2 pm emission. This suggests it should arise from a 
similarly balanced mixture of cations and neutrals from 
which the 17.8 pm band originates. The 17.4 pm band is 
predicted by the database to have similar emission inten¬ 
sities between cations and neutral FAHs, which is incon¬ 
sistent with our assignment solely to cations. Lastly, the 
database 15.4 pm band is seen to be cat ion-dominant, 
at odds with the neutral character of the 15.8 pm astro¬ 
nomical band from the observational data. It may sim¬ 
ply be the case that these two bands do not correspond 
to the same emission feature. Other possibilities for 
the observed inconsistencies may include variable PAH 
population abundances, incompleteness and/or poor rep¬ 
resentation of these particular (sub-)populations in the 
database, or simply that our data are not sufficiently sen¬ 
sitive to probe variations at the scale of two-to-one. An 
examination of the radial cuts shows the latter point to 
be true. 


5 . 1 . 2 . Molecular structure 

Molecular structure determines which vibrational 
modes are available to any particular PAH, but it is not 
known a priori how much of an influence structure will 
play in driving band intensity correlations on the whole. 


Some headway has been made in this area, in so far as 
particular band attributions have been suggested. The 

16.4 pm band has been attributed to FAHs containing 
pendent rings or “pointed” edges (e.g.. Van Kerckhoven 
et al. ]|2QQ01 [Peeters et al.||2QQ4| |Boersma et al.|| 2 U 10 f 
Ricca et al.|2010|). The former are single cycles that at- 
tach to the carbon skeleton through two shared carbon 
atoms, while the latter are corners that naturally arise 
in parallelogram-shaped FAHs, attached through three 
shared carbon atoms. The 16.4 pm band correlates well 
with the 17.8 pm band, which appears to be somewhat 
structure insensitive: the 17.8 pm ba nd was present in al l 
eight large neutral FAHs studied by Ricca et al. (2010), 
despite encompassing a variety of structures, including 
some with pointed edges. The 12.7 pm band, which cor¬ 
relates well with both the 16.4 and 17.8 pm emission, 
is associated with a combination of duo and trio C-H 


out-of-plan e modes (Hony et al.|2001 Bauschlicher et al 
2009 2010). Recently, this band was suggested to specit- 
ically originate in FAHs having an armchair-like struc¬ 
ture (jCandian et al.|2014|). As pendent rings and pointed 
edges naturally constitute duos and trios, it is then pos¬ 
sible that the 12.7 and 16.4 pm bands naturally correlate 
due to these common features. The 17.8 pm band is not 
known to have a dependence on duos or trios, but it has 
been suggested tha t it is a molecule-independent feature 
(|Ricca et al.||2010|), which may explain why it correlates 
with the 12.7 and 16.4 pm bands. 

The 17.4 pm band, which correlates with the 12.7, 

16.4 and 17.8 pm bands, generally appears to show 
the strongest emission intensities in com pact FAHs 
(Bauschlicher et al.||2010| Ricca et al. 2010). The 15.6 
pm band (which may or may not correspond to the astro¬ 
nomical 15.8 pm band) has been suggested to aris e from 
vibrations of pendent rings (Boersma et al.||2010). Pen¬ 
dent rings are also suggested as one possible origin for the 

16.4 pm band. If both the 15.8 and 16.4 pm bands are (at 
least in part) due to pendent rings, one would expect an 
intensity correlation, depending on environmental condi¬ 
tions (which we cannot speculate on here). In general 
the observational data show no such correlation between 
these bands, though there is a slight hint of a linear cor¬ 
relation when the ratio of 16.4/11.2 exceeds 0.12 (see Fig¬ 
ure^. The correlation of the 15.8 pm band with the 11.2 
pm oand is suggestive of large molecules, but it is known 
that the emission intensity of the 15.8 pm band decreases 
with increasing PAH siz e in the NASA Ames PAH IR 


Spectroscop i c Database dBauschlich er et al.||2QlQ| 
e£^y 201^ Boersma et^ al. ' 2014 bp . 1 ^'inally, |b c 

■. PUM ^ ^ ^ 


Ricca 


3 oersma 


et al.[ (| 20 10 |) examined classes of FAHs, including com- 
pact, irregular, and pent agon-containing FAHs. Minor 
systematic trends were observed, including the fact that 
irregular FAHs tend to have richer and more variable 
spectral features when compared to compact FAHs. As 
molecule size increased, compact FAHs exhibited little 
variation in the types of C-C-C vibrational modes that 
were active, whereas the complex edge structure of ir¬ 
regular FAHs fac ilitated a larger vari ety of C-C-C vi¬ 
brational modes (Boersma et al. 2010). Even with this 
knowledge, at the present there is no clear mapping be¬ 
tween molecular structure and the observational data. 


5 . 1 . 3 . Summary 
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The exact balance between charge and molecular struc¬ 
ture in driving the observed correlations of the 15-20 
|Lim bands cannot be discerned. However, the maps and 
radial cuts appear to cleanly separate these bands into 
those dominated by cations (11.0, 17.4 ttm), neutrals 
(11.2, 15.8 pm, 15-18 pm plateau), or a combination of 
the two (e.g., 12.7, 16.4, 17.8 pm bands). This suggests 
charge is the dominant factor driving these correlations, 
though there are indications that some bands should cor¬ 
relate due to structural arguments (e.g., the 12.7 and 16.4 
pm bands). Ultimately, theory points to wards a mixture 
of hig hly variable emission features (e.g., Boersma et al. 
2010| their Figures 11-15), whereas in space we observe 
a relatively simple picture of four bands and one under¬ 
lying plateau. This simplicity may reflect a basic com¬ 
monality of the carriers, or suggest that a relatively small 
collection of stable mol ecules are responsible for the ob¬ 
served emission bands (Boersma et al.||2 QlQ|). The latter 
is kn own as the gmndPAH hypothesis (Amdrews et al. 


5.2. Statistics of the 15-20 pm emission 

First we determine if the relative emission strengths of 
the 5-15 pm bands in our sample are consistent with sim - 
ilar measurements in the literature ( [Peeters et al. l2TOl) . 
We accomplish this by measuring the fractional contri- 
bution of the 6.2, 7.7, 8.6 and 11.2 pm bands to the total 
MIR emission (which we define here as the total flux 
emitted by these bands) for n ine Hll regions and fi fteen 
galaxies in our sample. Since Peeters et al. ([2002) only 
had two galaxies in their sample, we use the Hll regions 
as a direct point of comparison. Our results are presented 
in Figure 14, panel (a). We find that the H ll regions are 
consistent with those of [Peeters et al.|(|2002|), whose data 
is presented in panel (o), and there is no discrepancy 
bet ween Hll regions and galaxies in our sample. Note 
that Peeters et al. (2002[) included the 3.3 pm band and 
not the 8.6 pm band in their analysis, whereas we do the 
opposite (as our spectra start at 5 pm). The 3.3 and 
8.6 pm bands are quite weak (less than 10% of the total 
MIR emission in both studies) and should therefore not 
strongly influence the analysis. 

The same analysis is performed for the “FIR” fea¬ 
tures, i.e., the 15.8, 16.4, 17.4, 17.8 pm bands and the 
15-18 pm plateau, shown in Figure [l4{ panel (c). Due 
to the inconsistent detections of the ^.8, 17.4 and 17.8 
pm bands, the 16.4 pm band and the 15-18 pm plateau 
are used together as a proxy for the total FIR PAH 
flux. The 15.8, 16.4 and 17.8 pm bands have frac¬ 
tional strengths in Hll regions of 7 ± 4%, 16 ± 2% and 
6±1%, respectively. In galaxies, these are systematically 
lower, with values 4 ± 1%, 11 ± 3% and 3 ± 1%, respec¬ 
tively. The 17.4 pm band is observed to have similar 
strengths between H ll regions and galaxies (3 ± 1% in 
each). The plateau makes up the remainder of the mea¬ 
sured fractional strengths. The observed inconsistency 
by object type may be due to underestimating the 17.4 
pm PAH emission in Hll regions or overestimating the 
17.4 pm emission in galaxies. Neither would be surpris¬ 
ing, as Ceo emission at 17.4 pm is common, and it can 
only be removed by measuring the 18.9 pm Ceo emission, 
which itself is blended with the [S ill] 18.71 pm line. Re¬ 
call the typical 17.4/ 18.9 pm Can band strength ratio 
is approximately 0.5 ( [Garni et al.[[2QlQ Sellgren et al. 


2010 Bernard-Salas et al. 2012). Another possibility is 
that the 15.8, 16.4 and 1 y. 8 pni bands are all overesti¬ 
mated in H II regions, which requires a trio of systematic 
errors that we deem unlikely, as the same method was 
applied to both H ll regions and galaxies. The third pos¬ 
sibility is that these bands truly are systematically higher 
in Hll regions than galaxies. 

As a separate analysis, we investigate whether the 
MIR/FIR ratio depends on object type, where the 11.2 
and 16.4 pm bands are used as proxies for the MIR 
and FIR emission, respectively. We therefore examined 
490 low-resolution spectra of the SAGE-Spec sample and 
used the following selection criteria: 1. PAHs must be 
clearly present; 2. the spectra displayed no discontinu¬ 
ous jumps; 3. sources were less extended than 9"(which 
corresponds to the approximate size of the SL aperture 
at 6.2 pm); and 4. the 11.2 and 16.4 pm PAH bands 
must have 3cr detections. We identify 10 Hll regions and 
22 evolved stars as meeting these criteria. The Hll re¬ 
gions have a mean 16.4/11.2 ratio of 0.09 ±0.02. Evolved 
stars exhibited a ratio of 0.12 ± 0.05, consistent with the 
Hll regions. As such, we find no dependence on object 
type for the relative MIR and EIR band intensities in 
this analysis. 


6. SUMMARY & CONCLUSIONS 

The 15-20 pm PAH emission bands have been inves¬ 
tigated in fifty-seven sources, comprised of EMC point 
sources from the SAGE-Spectroscopy survey, nearby 
galaxies from SINGS, two extended cirrus sources and 
a spec tral map of NGG 7023 . We also included the sam¬ 
ple of [Boersma et al.[ ([2010[) and spectral maps of NGG 
2023 (papers 1 and 11). We investigated correlation plots 
of band flux ratios and examined the map morphologies 
of the RNe NGG 7023 and NGG 2023. We performed ra¬ 
dial cuts across these maps to evaluate the spatial profiles 
of the 15-20 pm PAH emission. Our primary conclusions 
are as follows: 


1 . Gorrelation plots: 

• The 16.4, 17.4 and the 17.8 pm bands are 
inter-correlated in all environments. 

• Within the RNe we see that the 16.4, 17.4 and 
17.8 pm bands further correlate with the 11.0 
and 12.7 pm bands. In Hll regions, the 12.7 
pm band also correlates with the 16.4 and 17.8 
pm bands. 

• The 11.2 and 15.8 pm bands and the 15-18 
pm plateau are inter-correlated in NGG 7023, 
consistent with results for NGG 2023 (paper 
I). The 15.8 pm band does not correlate with 
the 11.0, 12.7, 16.4, 17.4 or 17.8 pm bands 
in general. Within a sub-region of NGG 7023 
the 15.8 pm band appears to show a modest 
correlation with the 16.4 and 17.8 pm bands. 

2 . Spectral maps: 

• The maps show similar morphologies between 
the 11.0 and 17.4 pm bands, the 12.7, 16.4 and 
17.8 pm bands, and the 11.2, 15.8 pm bands 
with the 15-18 pm plateau. 
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Figure 14. Statistical summary of the 5-15 am and 15-20 
am PAH bands. The fraction of total flux in each band is illus¬ 
trated by its mean (white diamond for galaxies, black square for 
Hll regions), standard deviation (black line), and minimum and 
maximum values (denoted by the grey rectangles), (a) The 6.2, 
7.7, 8.6 and 11.2 a™ emission (called “MIR” here) are normalized 
to their combined total emission, (b) The same analys is, for class 
A and B sources, as presented in|Peeters et al.| (|2002|). Note the 
authors included 3.3 a^^^ emission, but did not have 8.6 a^^^ emis¬ 
sion. (c) The fractional emission in the “FIR” bands (15.8, 16.4, 
17.4, 17.8 aiii bands) are compared. The total FIR emission is con¬ 
sidered to be the sum of the 16.4 \im and plateau emission. The 
plateau fraction (> 0.7) is omitted for clarity. 


• The radial slices in the RNe maps show that 
the spatial profiles of the 12.7, 16.4 and 17.8 
pm bands can be reconstructed by averaging 
the 11.0 and the 11.2 pm band profiles. 

3. Charge dominates: we attribute the 12.7, 16.4 and 
17.8 pm bands to a combination of neutral and 
cationic molecules. The 17.4 pm band we attribute 
to cations, based on similarity of its spatial distri¬ 
bution to the 11.0 pm emission in the maps and 
radial cuts. Our results support the association of 
the 15.8 pm band with PAHs, specifically neutral 
PAHs (paper I). 

4. The 16.4/11.2 ratio in the general SAGE-Spec sam¬ 
ple shows no dependence on object type. The frac¬ 
tional emission of the 15.8, 16.4 and 17.8 pm bands 
in H II regions are systematically higher than those 
in galaxies. The 17.4 pm band is the exception, 
which has similar fractional strengths in both en¬ 
vironments. 


With the significant number of sources in our sample, 
we were able to discover correlations between the 15-20 
pm features that were not discernible in t he much smaller 
sample studied by Boersma et al. (2010). The forthcom¬ 
ing JWST instrurnent IVllKl will be an essential tool in 
furthering this analysis. Its large collecting area and in¬ 
creased spatial sensitivity will permit the construction of 
larger sample sizes, permitting further statistical analy¬ 
sis. Additionally, its improved spectral sensitivity will 
also allow us to examine the PAH band profiles in detail, 
with less line-blending than is possible with current data 
sets. These key improvements, in conjunction with theo¬ 
retical and laboratory analysis, may help disentangle the 
inter-related behaviours of the 15-20 pm bands. 


ACKNOLWEDGEMENTS 

We thank the anonymous referee for the useful feed¬ 
back. The authors thank J.D. Smith for supplying 
the SINGS data, as well as G. Boersma for supply¬ 
ing the data of their study. The authors also thank 
P.M. Woods for providing object classifications. The au¬ 
thors acknowledge support from NSERG discovery grant, 
NSERG acceleration grant and ERA. The IRS was a col¬ 
laborative venture between Cornell University and Ball 
Aerospace Corporation funded by NASA through the 
Jet Propulsion Labo ratory and Ames Research Center 
(Houck et al. 2004). This research has made use of 
NASA's Astrophysics Data System Bibliographic Ser¬ 
vices, and the SIMBAD database, operated at CDS, 
Strasbourg, Erance. 


REFERENCES 


Allamandola, L. J., Hudgins, D. M., & Sandford, S. A. 1999, ApJ, 
511, L115 

Allamandola, L. J., Tielens, A. G. G. M., Sz Barker, J. R. 1989, 
ApJS, 71, 733 

Andrews, H., Boersma, C., Werner, M. W., et al. 2015, ApJ, 807, 
99 

Baschek, B., Beltrametti, M., Koppen, J., & Traving, G. 1982, 
A&A, 105, 300 

Bauschlicher, Jr., C. W., Peeters, E., & Allamandola, L. J. 2009, 
ApJ, 697, 311 
































































Probing charge with the 15-20 jam PAH bands 


17 


Bauschlicher, Jr., C. W., Boersma, C., Ricca, A., et al. 2010, 
ApJS, 189, 341 

Beasley, A. J., Ellingsen, S. P., Claussen, M. J., & Wilcots, E. 
1996, ApJ, 459, 600 

Bernard-Salas, J., Cami, J., Peelers, E., et al. 2012, ApJ, 757, 41 
Bernard-Salas, J., Peelers, E., Sloan, G. C., et al. 2009, ApJ, 699, 
1541 

Berne, O., & Tielens, A. G. G. M. 2012, PNAS, 109, 401 
Boersma, G., Bauschlicher, G. W., Allamandola, L. J., et al. 2010, 
A&A, 511, A32 

Boersma, G., Bauschlicher, Jr., G. W., Ricca, A., et al. 2011, 

ApJ, 729, 64 

Boersma, G., Bregman, J., & Allamandola, L. J. 2014a, ApJ, 795, 
110 

Boersma, G., Bregman, J. D., &: Allamandola, L. J. 2013, ApJ, 
769, 117 

Boersma, G., Bauschlicher, Jr., G. W., Ricca, A., et al. 2014b, 
ApJS, 211, 8 
—. 2015, ApJ, 806, 121 

Boulanger, F., Boisssel, P., Gesarsky, D., &: Ryter, G. 1998, A&A, 
339, 194 

Garni, J., Bernard-Salas, J., Peeters, E., & Malek, S. E. 2010, Sci, 
329, 1180 

Gandian, A., Sarre, P. J., Sz Tielens, A. G. G. M. 2014, ApJ, 791, 
LIO 

Devlin, S. J., Gnanadesikan, R., & Kettenring, J. R. 1975, 
Biometrika, 62, 531 
Finkenzeller, U. 1985, A&A, 151, 340 

Galliano, F., Madden, S. G., Tielens, A. G. G. M., Peeters, E., Sz 
Jones, A. P. 2008, ApJ, 679, 310 
Guetter, H. H. 1968, PASP, 80, 197 
Habing, H. J. 1968, BAN, 19, 421 

Hony, S., Van Kerckhoven, G., Peeters, E., et al. 2001, A&A, 370, 
1030 

Houck, J. R., Roellig, T. L., van Gleve, J., et al. 2004, ApJS, 154, 
18 

Hudgins, D. M., Sz Allamandola, L. J. 1999, ApJ, 516, L41 
Kemper, F., Woods, P. M., Antoniou, V., et al. 2010, PASP, 122, 
683 


Kennicutt, Jr., R. G., Armus, L., Bendo, G., et al. 2003, PASP, 
115, 928 

Mookerjea, B., Sandell, G., Jarrett, T. H., Sz McMullin, J. P. 
2009, A&A, 507, 1485 

Moutou, G., Verstraete, L., Leger, A., Sellgren, K., Sz Schmidt, 

W. 2000, A&A, 354, L17 

Peeters, E., Hony, S., Van Kerckhoven, G., et al. 2002, A&A, 390, 
1089 

Peeters, E., Mattioda, A. L., Hudgins, D. M., & Allamandola, 

L. J. 2004, ApJ, 617, L65 

Peeters, E., Tielens, A. G. G. M., Allamandola, L. J., Sz Wolfire, 

M. G. 2012, ApJ, 747, 44 

Ricca, A., Bauschlicher, Jr., G. W., Mattioda, A. L., Boersma, G., 
& Allamandola, L. J. 2010, ApJ, 709, 42 
Rosenberg, M. J. F., Berne, O., Boersma, G., Allamandola, L. J., 
& Tielens, A. G. G. M. 2011, A&A, 532, A128 
Sellgren, K., Uchida, K. I., & Werner, M. W. 2007, ApJ, 659, 1338 
Sellgren, K., Werner, M. W., Ingalls, J. G., et al. 2010, ApJ, 722, 
L54 

Sheffer, Y., Wolfire, M. G., Hollenbach, D. J., Kaufman, M. J., Sz 
Gordier, M. 2011, ApJ, 741, 45 

Smith, J. D. T., Armus, L., Dale, D. A., et al. 2007a, PASP, 119, 
1133 

Smith, J. D. T., Draine, B. T., Dale, D. A., et al. 2007b, ApJ, 
656, 770 

Stock, D. J., Peeters, E., Ghoi, W. D.-Y., & Shannon, M. J. 2014, 
ApJ, 791, 99 

Stock, D. J., Peeters, E., Tielens, A. G. G. M., Otaguro, J. N., & 
Bik, A. 2013, ApJ, 771, 72 

Tielens, A. G. G. M. 2005, The Physics and Ghemistry of the 
Interstellar Medium 
—. 2008, ARA&A, 46, 289 
—. 2013, RvMP, 85, 1021 
Treffers, R., Sz Gohen, M. 1974, ApJ, 188, 545 
Van Kerckhoven, G. 2002, PhD thesis. Institute of Astronomy, 
Katholieke Universiteit Leuven, Belgium 
Van Kerckhoven, G., Hony, S., Peeters, E., et al. 2000, A&A, 357, 
1013 

Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 
1 



